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Abstract : 

T— I ■ 

This work is devoted to the study of a particular four-body quantum problem of Wolfes type on the 
line [1]. We solve the corresponding Schrddinger equation in an harmonic confinement for the particles 
where we provide the full eigensolutions and corresponding eigenenergies. We discuss also the domain 
of the coupling constant where the irregular solution becomes physically acceptable. Finally, we extend 
the problem to D-dimensional space. 



PACS: 02.30.Hq, 03.65.-w, 03.65.Ge 



1 Introduction 



There exists a very limited number of exactly solvable many-body systems, even in one dimension 
(ID) [2, 3]. The Calogero model constitutes one of the famous ones which was exhaustively studied 
[4, 5]. A survey of many quantum integrable systems was done by Olshanetsky and Perelomov [6]. 
They classified the systems with respect to Lie algebras. Point interactions have also been considered, 
still in D = 1 dimensional space [7, 8]. 

The quest for exactly solvable non trivial quantum problems of few interacting particles on the line still 
retains attention. Early works of three-body linear problems of Calogero-Marchioro-Wolfes [9, 10, 11] 
have been followed by new extensions and cases. In a non exhaustive way, we can quote the three-body 
problems of Khare et al. [12], Quesne [13] and Meljanac et al. [14]. 

In a recent paper, we have exactly solved a generalization of the three-body Calogero-Marchioro- 
Wolfes problem with an additional non-translationally invariant three-body potential [15]. In this 
paper we propose a "completely solvable" model, i.e., a model were we give explicitly the eigenvalues 
and the complete set of eigensolutions of the problem. This completely solvable four-body bound 
states problem, in one and D-dimensional space, can be considered as a generalization of the four- 
body scattering linear problem proposed and solved by Wolfes [1] with additional non-translationally 
invariant four-body potentials. 

The paper is organized as follows. In section 2 we expose and solve the problem for the harmonic case. 
The section 3 is devoted to extension to D-dimensional case. Our conclusions are drawn in section 4. 



2 A generalization of the linear Wolfes four-body problem 



We consider the four-body Hamiltonian on the line 
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This Hamiltonian represents a system of four particles on the line with the same mass (with units 
h = 2m = 1) interacting via only four-body potentials. One potential is translationally invariant 
with coupling constant A of Wolfes type [1] and the two others with coupling constants fi, j3 are not 
translationally invariant. The whole system is confined in an harmonic oscillator trap. 



The problem is solved in the following way. Setting 
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the Hamiltonian, Eq.(2), becomes 
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Note that, if /3 = 0, the problem is separable and the derivation of the solutions is straightforward. 
Otherwise, it is not separable in {R, s, t, u} variables. To overcome this situation we introduce the 
following hyperspherical coordinates: 



R = rcosa, s = r sin a cos t = r sinasmO sincp, u = r sin a sin 9 cos ip, 
< r<oo, 0<a<vr, < 6» < vr, < < 2vr. 



(5) 



The stationary Schrodinger equation is then written as: 
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V(r,a,0,<p) =EV(r,a,9,(p), 



where *(r, a, 9, (p) represent the eigensolutions associated to eigenenergy E. 



(6) 



This four-body problem described by this equation (6) may be mapped to the problem of one particle 
in four dimensional space, with a non central potential of the form 



V(r,a,6,<p) = f 1 (r) + \(f 2 (a) + 
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It is then clear that the problem becomes separable in the four variables {r,a,9,<p}. To find the 
solution we factorize the wave function as follows: 
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Accordingly, equation (6) separates in four decoupled differential equations: 

\--r~2 ^ — ^2 — ~ ^ *~T~ I = D £mn Ge mn (a), (11) 
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In the interval < <p < 2ir the potential involved in equation (9) has a periodicity of tt and has 
singularities at ip = kir/2,k = 0,1,2,3. The equation is first solved in the interval [0, 7r/2]. It can 
be treated if and only if A > —1/4 otherwise the operator has several self-adjoint extensions, each of 
them leading to a different spectrum [16, 17]. The regular eigensolutions of equation (9) on [0, tt/2] 
read [6, 18] 

<& n (^) = (sin2 9? ) 1 / 2+a C^ 2+a \cos2 V ) (13) 
where the Cn^ 2+a ^ denote the Gegenbauer Polynomials [19]. The corresponding eigenvalues are 

B n = 4 ( I + a + n) (14) 
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with 



n = 0,1,2,... a= Vi+iA (l5) 

The extension to the whole interval [0, 2tt] is made in two steps. First one, from [0, tt] up to [tt, 2tt] 
using the periodicity of the solutions. On the other hand, the symmetric extension of the solutions in 
the interval [tt/2, tt] reads 

$ n (<p) = (- sin2 9 9) 1 / 2+a Cl 1/2+a) (cos29?) (16) 
so that the real power of (— sin2y?) is defined. In a compact form we obtain 

if> n (<p) = (e lS in2^) 1 / 2+a C^ 2+a \cos2<p) (17) 
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Note that when a = 1/2, (A = 0), a 5-pathology occurs at if = tt/2. 

On the other hand, the antisymmetric extension of the solution reads 

$ n (<p) = sgn(sin(2 V 9))(eisin2o 5 ) 1 / 2+a C T ( t 1/2+a) (cos 2^) (18) 
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where sgn(x) = x/\x\,x / 0, denotes the sign of the variable x. 



Generally, only the regular solution, $n + \ corresponding to 1/2 + a, is retained. However, the irregular 
solution, § n \ which is distinct from <E>n + ^ for A > —1/4, corresponding to 1/2 — a is physically 
acceptable when the Dirichlet condition is satisfied for —1/4 < A < (attractive potentials). If we 
release the Dirichlet condition, and ask only for the square integrability of the solution, as in [20], then 
$1 ^ can be retained for —1/4 < A < 3/4 [15]. For the irregular antisymmetric solution, a (derivative 
of) ^-pathology occurs at if = ir/2, when a = 1/2, (A = 0). 

Introducing the parameter 

b n = J~B~ n = 1 + 2o + 2n (19) 

(where we have only considered the positive root b n = \/B n because the other root b n = —y/B n leads 
to non-square integrable solutions for most values of n) the solution of Eq.(10) on the whole interval 
[0, 7r] can be obtained as [6, 15] 

@m,n( 9 ) = sgn(cos60 se (sin#) 6 " +1 / 2 (e 2 cos#) c+1/2 P^™' c) (cos 20) (20) 
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corresponding to the eigenvalue 

C m . a = {2m + b n + c+lf . (22) 

In fact, a = c, if one takes the definitions Eqs.(15,21). In Eq.(20), the Pm n ' c ^ denote the Jacobi 
polynomials [19]. The value s$ = 0, 1 according to the fact that the solution has been extended in a 
symmetric (antisymmetric) way from 6 € [0, tt/2] to 6 G [tt/2,tt]. Note that, when c = 1/2, (A = 0), a 
5-pathology occurs for the symmetric solution at 6 = tt/2. 

The equation (11) is solved in the same manner as for equation (10). Setting 



Cm,n = \fC m>n =2m + b n + c+l (23) 



(where we have only considered the positive root c m ^ n = yJC m , n because the other root c m , n = —^/C„ 
leads to non-integrable solutions for most values of m) the solution writes as 



Gt,m, n (a) = sgn(cosa) Sa (sina) Cm '™ +1/2 (e 3 cos a) d+1/2 m ' n ' d) (cos2a) (24) 
£ = 0,1,2,... d=iyiT4^ (25) 

, -, 1 - e 3 7T 3 - e 3 vr 

^3 = ±1 — ^ — — < a < — - — — 
2 2 ~ ~ 2 2 

corresponding to the eigenvalue 

D e ,m, n = (2^ + c m , n + d + l) 2 . (26) 

The value s a = 0, 1 according to the fact that the solution has been extended in a symmetric (anti- 
symmetric) way from a £ [0,7r/2] to a G [tt/2,tt]. Note that when d = 1/2, (/x = 0), for symmetric 
solutions, a (5-pathology occurs at a = 7r/2. 



On the other hand, the reduced radial equation (12) is solved in the interval < r < oo, with the 
condition of square integrability for the solutions. It implies i^, m , n (r) — > as r — > oo. 

We have to impose /3 + Di^m,n > in order to treat the centrifugal barrier in the vicinity of r — 0. 
Note that taking (3 + D^ m n = leads to several self-adjoint extensions parametrized by a phase. 
This fact has been discussed in [15]. More details can be found in [21, 22]. Also attractive barriers 
/3 + £V,m,n < 0, mentioned in [15], have been treated in [23, 24, 25]. Taking into account the definition 
of -D^ min , Eq.(26), we must have 

P + D t ^ n = (3 + (2l + 2m + 2n + 2a + c + d + ?>) 2 > 0, Vn > 0, Vm > 0, W > 0. (27) 

For positive values a, c, d, the quantity /3 + -D^ m ,n is minimal for n = 0, m = 0, 1 = and a = c = d = 
( recall that a > 0, see (15), that c > 0, see (20) and that d > see (25)). It puts constraint on j3 to 
satisfy (3 > —9 when a,c and d are positive. We introduce the auxiliary parameter Kt^ m ^ n defined by 

K ?,m,n = P + A,m,n, «*,m,n = ^//? + A?,m,n • (28) 

The solution of the radial equation (12) is [15] 

F M , m , n (r) = r^,n+± exp (-^pj L|f V 2 ), fc = 0, 1, 2..., (29) 

and it is associated to the eigenenergy 

£ M , m , n = 2w(2A; + Ke, m , n + 1), A: = 0, 1,2.... (30) 
The are the generalized Laguerre polynomials [19]. 



Taking into account all information, we conclude that the physically acceptable solutions of the 
Schrodinger equation (6) are 
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Here s 2(/p = 0, 1 according to the parity of the solution for 92 G [0, 7r]. We note that ^-pathologies occur 
for respectively d = 1/2, a = c = 1/2, in Eq.31, when symmetric solutions ( i.e., s a , sg or s 2((3 equal to 
zero) are considered. 



The normalization constants Nk £ m n can be calculated from 
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Use is made, here, of the orthogonality properties of Gegenbauer, Jacobi and Laguerre polynomials 
[26]. 

The full expression of the eigenenergies is expressed by 

E k/ , m ,n = E kt2e+2m+2n = 2Lo(2k + l + ^ + (2£ + 2m + 2n + 2a + c + d + 3) 2 ) , (33) 
k = 0,1,2,..., £ = 0,1,2,... m = 0,l,2,... n = 0,l,2,..., . 

For illustration, the equation (31), multiplied by r 3 / 2 (sina)V sin 9 reads in Cartesian coordinates, for 
the symmetric case (we remind the reader that a = c = + A) 

Vl/(xi,X 2 ,X3,X 4 ) (X |(xi + X 4 - X 2 - X 3 )(xi + X 3 - X 2 - X 4 )(xi + X 2 - X 3 - X 4 )| a+1 / 2 

„ [(*, - I2 f + (x 3 - X 4 )T C»«/ 2 ( 2(»! - *,)(*.-*») 

V(xi - x 2 r + (^3 - Xa) Z , 

x [4(x 2 + x 2 . + xl + xl) - (xi + x 2 + x 3 + x 4 ) 2 ] m 

x p(2n+2o+l,o) f 8(XjX 2 + X 3 X 4 ) - (x 1 + X 2 + X 3 + X 4 ) 

m l^4(xf + x\ + xl + x\) - (xi +x 2 + x 3 + x 4 ) 2 y 

x |xi + x 2 + x 3 + x 4 | d+1/2 (x 2 + x 2 2 + x^ + xl) 9 

p (2m+2n+3a+2,d) ( 2 ( x l + ^2 + ^3 + x a) 2 - (xf + xj + xj + x|) \ 

* \ xf + x 2 . + xfj + x 2 y 

x exp(-w(x 2 + x^ + x^ + xl)/2) L^(w(x 2 + x^ + x^ + xl)) (34) 

p = + (21 + 2m + 2n + 3a + d + 3) 2 
q = p /2- (2m + 2n + 3a + d + 3)/2 . 

The eigenvalue Eq.(33) shows a degeneracy, i.e., all solutions such that the equality I + m + n = N 
is satisfied, iV being fixed, correspond to the same eigenvalue. This imply that any combination of 
solutions such that l + m + n is constant is also solution. We can then obtain, k,£,m + n being fixed, 
a unique solution, symmetric in the permutation on the set of variables {x±, x 2 , x 3 , x 4 } when n + m is 
even. This latter symmetric solution is in fact proportional to J2<reS 4 ^( x c(i)i x a{2)-> x cj(3)-, x o(a)) where 
Si denotes the permutation group of four elements. 

Let us now consider the irregular "polynomial" solutions corresponding to 1/2 — a. We have to replace 
a by —a in all equations, from Eq.(17) until Eq.(33). Recall that for —1/4 < A < 3/4, the irregular 
solutions, Eq.(17), are square integrable, as seen before. The Sturm-Liouville operator (9) is self- 
adjoint for A / —1/4. It has to be added that, a S pathology occurs in (31) for a = 1/2 (A = 0). We 
consider values of A in ] - 1/4, 0[U]0, 3/4[. 



with 



We next examine the impact on the change a i-> —a on the function m] „(9), Eq.(20). We remind 
the reader that |c| = \a\,a = ±>/l/4 + A. Here a = — + A. Also, we allow c = ±y/l/4 + A 
i.e., we consider negative values of c. We have \a\ = \c\ but we allow a = ±c (both quantities 
being allowed to be negative and having a different sign). Note that another situation corresponds to 
a = + A = — c. But now we concentrate on a = —y/T/4 + A, c = ±a. 

First, to ensure the self-adjointness of the Sturm-Liouville operator Eq.(lO), we have to impose both 
B n 7^ for every n and A 7^ —1/4. Taking into account Eq.(19), for a changed into —a, we have 
clearly to ask for a 7^ 1/2 which is equivalent to A / 0. These latter conditions are fulfilled for A in 
]-l/4,0[U]0,3/4[. ' 

Then we examine for which value of a we obtain square integrable solutions G m , n (0) for every value 
of n. Since 

^ + b n = 2n + ^-2a (a = ^v / l + 4A) , (35) 

we have 

(Vn > 0) 1 + 6 n > I - 2a . (36) 

The function m , n (0), Eq.(20), leads to square integrable solutions for every value of n if a < 1, 
corresponding to b n + 1 > 0, (Vn). This latter inequality happens for A < 3/4. This is satisfied when 
A e] — 1/4, 0[U]0, 3/4[ (see the paragraph below Eq.(17)). As far as the term A/ cos 2 is concerned 
irregular solutions happen when c = ^/l/4 + A is changed in — c in Eq.(20). These solutions are square 
integrable when c < 1 i.e., for A < 3/4. 

Consider the last angular equation concerning the variable a. The differential operator is self-adjoint 
provided that 

(Vm>0)(Vn>0) c nhTl = 2m + 2n + 2 ± c - 2a > 2 - 2a ± c , (37) 

is non zero. So that the self-adjointness is ensured when 2 — 2a ± c > or equivalently 2 — 2a ± a > 0. 
The quantity 2 — a is always positive in the domain of acceptable A whereas 2 — 3a is positive for 
A < 7/36 ~ 0.194444. Also a has to be different from -1/4. 

The square integrability of the function G^ )m , jra , Eq.(24), is ensured by c m) „ + 1 > 0, (Vm) (Vn) i.e. 
when 3 — 2a ± c > or equivalently 3 — 2a ± a > . Both quantities 3 — 2a ± a are alway positive 
for —1/4 < A < 3/4. This defines a domain in A of acceptable solutions. Also if d is allowed to be 
negative we must have — 1/4 < d < 3/4. 

As far as the radial equation is concerned, the constraint j3 + D^ m ^ n > allows us to treat the 
centrifugal barrier in the vicinity of r = 0. (see the above discussion for j3 + T>^ mn < 0.) Taking into 
account the definition of Z^, m , n , Eq.(26), we have 

/3 + A?,m,n = /3 + (2£ + 2m + 2n-2a±c + d + 3) 2 > 0, Vn > 0, Vm > 0, W > 0. (38) 

This is satisfied for every {£, m, n} when (3 > and (see Eq.(38)) 

±c-2a + 3 + d> /3<0 (39) 

This condition defines a domain of acceptable values of (3 depending on the values of A reminding that 
A e] — j,0[L)]0, |[. Under such conditions, the radial solutions, Eq.(29), are square integrable because 

K£,m,n > 0. 



Note that the spectrum for irregular solutions has eigenvalues lower than the ones corresponding to 
the regular solutions. This spectrum is given by 



2oj (2k + 1 + y/p+(2£ + 2m + 2n±c + d-2a + 3) 2 ^) , (40) 
k = 0,1,2,..., £ = 0,1,2,... m = 0,l,2,... n = 0,l,2,..., . 
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Also d can be changed in — d = — a/1/4 + fi. The requirement 



±c - 2a + 3 - d > V^fi j3<0 

±c-2a + 3-d >0 /3>0 (41) 



ensures the self-adjointness of Eq.(12). 



3 D-dimension 



The generalization to the D-dimensional space follows the same strategy as in section 2. We consider 
the Hamiltonian: 



ff = E(-A j+ u,Vf) + 2A £ 1 + «f^ + (42) 
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Setting as in section 2 
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the Hamiltonian Eq.(42) becomes 



i? = -A^ — Ag — Af- A a + uj 2 {R 2 + s 2 + t 2 + u 2 ) 

+ 1$ + A (^2 + ^ + ^2) + #2 + s 2 + t 2 + u 2 • ( 44 ) 



Since the potential does not depend on the angles between i?, s, t, u, the separation of angular and 
radial variables, together with the use of the hyperspherical harmonics [19], allows us to write 



= (fl^C5=^ *WW n *) ^,w(^)^,[M f ](^)^,[M„](^) • (45) 

Here [M] denotes the set [M\ = {mi,m2, m p },p = D — 2 satisfying £ = my > mi > m,2~- > m p > 
and [19] 

^,[M] = ****** n(sm^r* ncir t/2 (»»w) (46) 

k=l fc=0 

where the hyper spherical harmonics Yi\m\ are given in terms of the Gegenbauer polynomials C%(x). 
For D = 2, Yi = exp(i£(j)),£ G Z and G [0, 2tt]. Recall that the hyperspherical polar coordinates are 
given by 



xi = r cos 9\ 

X2 = rsin^x cos# 2 

X3 = r sin 9\ sin O2 cos #3 

x p+ i = r sin 9\ sin 62 ■ ■ sin 6 P cos </> 

= rsin^i sin 62-. sin 6 P sin (ft (47) 

with 6> fc , G [0,7r], fc = 1,2, ..p and G [0,27r]. Setting (D - 3)/2 = md (here > -1/2) we obtain 



d 2 d 2 d 2 d 2 u+(£ R + md)(£ R + md+l) 



\ dR 2 ds 2 dt 2 du 2 K ' R 2 

A + (£ s + md)(£ s + md+l) X + (£ t + md)(£ t + md+l) 
2 ' 
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t 2 



+ A + (4 + m ^" + md+1) + ^ + / + ^ ) ^^(^.-.M) - (48) 

We introduce the following hyperspherical coordinates: 

R = r cos a, s = r sin a cos 6* , t = r sin a sin 9 sin ip, u = r sin a sin 6 cos 99, 
< r<oo, 0<a<|, < < |, < y> < |. (49) 

Here a,0,<p G [0,7r/2] because R,s,t,u are positive. 

The stationary Schrodinger equation is then written: 
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V(r,a,0,<p) =E^(r,a,0,ip), 



(50) 



where ty(r,a,8,(p) represents the eigensolution associated to eigenenergy E. 
Assuming 

*M,™M r > «, 0, if) = '->-}= '-^ ^=-^ n (ip) 51 

we obtain 

/ d 2 X + (i t + md)(i t + md + 1) X + (i u + md)(i u + md+\)\ 
V d<p z sm z (p cos 2 p I 
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and 
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aar sir Q COS 2 Q ' 



+ w2r2 H ^T^ 1 ) F k,£,m,n( r ) = E k,e,m,n F k/ ^ n (r) . (55) 



The equation (52) is solved in the interval [0,7r/2]. The regular eigensolutions of equation (52) on 
[0,tt/2] read [6, 18] 

= (sin^) 1/2+a (cos(^) 1/2+fc P,( a ' b) (cos2v9) (56) 
where the pjf'^ denote the Jacobi Polynomials [19]. The corresponding eigenvalues are 

B n = (a + b + 1 + 2n) 2 (57) 

with 

n = 0,1,2,... a = ^Jx+(i t + md+ 1/2) 2 b = ^X + (i u + md + 1/2) 2 (58) 

The operator in Eq.(52) is self-adjoint for every it,i u provided that A + (md + 1/2) 2 > 0. Introducing 
the parameter 

b n = ^/B^ =l + a + b + 2n (59) 

(where we have only considered the positive root b n = \[B~ n because the other root b n = — \/B^ leads 
to non-square integrable solutions for most values of n) the solution of Eq.(53) on the interval [0, tt/2] 



can be obtained as [6, 15] 

e ro ,„(0) = (sin0) 6 " +1 / 2 (cos#) c+1 / 2 p£ n > c \cos26) (60) 
m = 0,l,2,... c = + (4 + 1/2) 2 

corresponding to the eigenvalue 

C m , n = (2m + 6 n + c+l) 2 . (61) 

In Eq.(60) the Pm"'^ denote the Jacobi polynomials [19]. The operator in Eq.(53) is self-adjoint 
for every £ t ,£ u ,£ s provided that A + (md + 1/2) 2 > 0. The second condition namely B n ^ i.e. 
2^A + (md + 1/2) 2 + 1 > is always satisfied. 

The equation (54) is solved in the same manner as done for equation (53). Setting 



-m,n 



y/C m ,n =2m + b n + c+l (62) 



(where we have only considered the positive root c mi „ = y/C m>n because the other root c m . n = — ^/C m ,n 
leads to non-integrable solutions for most values of m) the solution writes as 

G t , m , n (a) = (sma) c ^ +1/2 (cosa) d+1/2 P £ (cm ' n,d) (cos 2a) (63) 



£ = 0,1,2,... d= yJv + (£ R + md+l/2) 2 (64) 

corresponding to the eigenvalue 

D e , m , n = (2£ + c m , n + d + l) 2 . (65) 

The operator in Eq.(54) is self-adjoint for every £t,£ u , £ S ,£ R provided that A + (md + 1/2) 2 > 0. The 
condition C m ^ n / i.e. 3\/A + (md + 1/2) 2 + 2 > is always satisfied. 

On the other hand, the reduced radial equation (55) is solved in the interval < r < oo, with the 
condition of square integrability for the solutions. It implies Fk/ ; m ;n (r) — > as r — > oo. We have to 
impose j3 + D^ m;n > in order to treat the centrifugal barrier in the vicinity of r = 0. 

P + A?,m,n = P + (21 + 2m + 2n + a + b + c + d + 3) 2 > 0, Vn > 0, Vm > 0, W > 0. (66) 

For positive a, b, c, d the quantity P + Dt mn is minimal for n = 0, m = 0, £ = and a = b = c = d = 
( recall that a, b > 0, see (58), that c > 0, see (60) and that d > see (64)) It put constraint on /3 to 
satisfy (3 > —9 when a, 6,c and d are positive. We introduce the auxiliary parameter Ke,m,n defined by 



„2 



(3 + D e ,m,n ■ 

(67) 

The solution of the radial equation (55) is [15] 

F kAm ,n(r) = r Ke ^ +l * exp 4^'" V'), k = 0,1, 2..., (68) 

and it is associated to the eigenenergy 

E kAm , n = 2u(2k + K t , m , n + 1), fe = 0, 1,2.... (69) 



The are the generalized Laguerre polynomials [19] . 

Taking into account all information, we conclude that the physically acceptable solutions of the 
Schrodinger equation (50) are 



X (Sin a )2 m +2n+a+6 +C +| (cog a) d + l/2 p (2m+2n+a+b+c+2,d) ^ ^ 
X (sin ef n+a+b+ i (COS 0y+^p^n+a+b+l,c) ( cog 

x (sin (^) a+ 5 (cos <f) h+h 2p( a ^ (cos 2<p), (70) 
k = 0,1,2,..., £ = 0,1,2,..., m = 0,1,2,...,, n = 0,l,2,..., 

The normalization constants N^^m^ can be calculated from 

r+oo /" 7r /2 /" 7r /2 

/ r 3 dr sin 2 (a) da / s'm 6 d6 dip ^k,e,m,n{r, a, 9, (p)^k',e',m',n'{r, a, 9, <p) 
Jo Jo Jo Jo 

= Sk,k'h,i'?>m,m'&n,riNk,l,m,n> (71) 

use is made, here, of the orthogonality properties of Gegenbauer, Jacobi and Laguerre polynomials 
[26]. 

The full expression of the eigenenergies is expressed by 

E k ,e, m ,n = E k/+m+n = 2w ^2fc + 1 + ^f3+(2£ + 2m + 2n + a + b + c + d + 3) 2 ) , (72) 
k = 0,1,2,..., £ = 0,1,2,... m = 0,1,2,... n = 0, 1, 2, . 



4 Conclusions 



In this paper, we have proposed and solved a linear four-body quantum problem of Wolfes type, in 
the case of bound states. The operator we consider is composed of an Harmonic trap with four-body 
translationally invariant potential proposed by Wolfes with additional non-translationally invariant 
four-body potential. We explicitly give the full solutions of the corresponding Schrodinger equation, 
namely the wave functions in terms of the angular and radial variables together with the energy spec- 
trum. 

We have been able, also, to generalize and solve this linear problem to any D-dimensional space where 
the explicit eigensolutions and the corresponding spectrum have been calculated. 
We also investigate the domain of coupling constant where the irregular solutions become square in- 
tegrable and physically acceptable. 

This proposed exactly solved four-body problem can be generalized to a "Coulomb-type" problem 
where the Harmonic trap can be replaced by an attractive "Coulomb" type potential as done in [12] 
for the three-body problem. 
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